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The distortion of a two-dimensional bubble (or drop) in a corner flow of an in viscid in¬ 
compressible fluid is considered. Numerical solutions are obtained by series truncation. The 
results confirm and extend previous calculations. 


1 Introduction 

Over the years many solutions have been obtained for the deformation of two dimen¬ 
sional bubbles in a flow. Both steady and unsteady solutions were considered. Elere we 
concentrate our attention to steady solutions, (see, for example, Baker & Moore [1], and 
the references cited in that paper for time dependent solutions). 

McLeod [4] obtained an exact solution for a free bubble rising in a fluid. Vanden- 
Broeck & Keller [12] showed numerically that McLeod’s solution is a member of a one 
parameter family of solutions. They also solved the flow past a bubble attached to a 
wall (i.e. the configuration of Figure 1 with a = ji). This problem was also considered by 
Meiron [5], Shankar [6] and Tanveer [7]. Vanden-Broeck & Keller [13] considered the 
related problem of the deformation of a bubble attached to the walls of a rectangular 
wedge (i.e. the configuration of Figure 1 with a = tc/ 2). When reflected into the walls, this 
configuration models the distortion of a bubble in a straining flow. 

The solutions in Vanden-Broeck & Keller [12, 13] were obtained numerically by 
boundary integral equation methods. It was found that the families of solutions ultimately 
approach limiting configurations in which opposite sides of the bubble touch each other. 
In Vanden-Broeck & Keller [13], numerical solutions were not obtained close to the 
limiting configuration, and an approximate analytical solution was derived to describe the 
limiting configuration. This approximate solution was extended to other values of a in 
Vanden-Broeck [8]. It was also shown [8] that solutions can be found past the limiting 
configurations by forcing contact of the opposite sides at one point or along a segment. 

In this paper, we present another numerical procedure to compute solutions for the 
flow configuration of Figure 1. The scheme is based on series truncation, and is similar 
in philosophy (if not in details) to those used by Vanden-Broeck [9-11] and Daboussy 
et al. [3]. Numerical solutions are obtained for all values of a and of the contact angle /(. 
The results of Vanden-Broeck & Keller [12, 13] are recovered for a = jt and a. = tt/ 2. In 
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Figure 1. Sketch of the flow and coordinates. 


addition, the results of Vanden-Broeck & Keller [13] are extended numerically up to the 
limiting configuration. 

The solutions for values of a not equal to n or n/2 are not covered in the calculations 
of Vanden-Broeck & Keller [12, 13], although approximate solutions based on slender 
body theory were derived [8]. The fully nonlinear solutions derived in this paper are of 
interest, since flows in wedges have recently been used to model the flow in a crack in a 
seawall struck by a wave [2]. The present calculations provide an additional model for 
cracks by allowing a trapped bubble at the apex of the wedge. 

The problem is formulated in § 2. The numerical procedure is described in § 3 and the 
numerical results are presented in § 4. 


2 Formulation of the problem 

We consider the two-dimensional flow configuration shown in Figure 1. The flow domain 
is bounded by the two walls AE and CD and by the free surface AC. The fluid is assumed 
to be inviscid and incompressible, and the flow to be irrotational and steady. We take 
into account the effect of the surface tension T at the interface, but we ignore the flow 
inside the bubble OACO, assuming that the pressure is a constant Pg throughout it. The 
angle between the two walls is denoted by a, and the contact angle between walls and the 
free surface is denoted by /? (see Figure 1). We introduce Cartesian coordinates with the 
origin at the intersection of the two walls and the x-axis along the wall CD. 

We introduce the potential function </> and the streamfunction ip. We choose tp = 0 
along the walls and along the bubble and r/j = 0 at the point B in the middle of the free 
surface profile (see Figure 1). We denote by +b the values of cj) at C and A. Next, we 
define dimensionless variables by taking T/(pb) as the reference velocity and pb 2 /T as 
the unit length. Then the flow region corresponds to the upper half (i.e. i p ^ 0) of the 
(0, i p) plane. The bubble surface corresponds to the segment — 1 < cp < 1, xp = 0 of the 
0 axis (see Figure 2). 
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Figure 2. The complex potential plane. 


On the free surface AC, Bernoulli’s equation yields 


i^+e-) T K = Ps - fB , 

P P 

(1) 

where u and v denote the x and y components of the velocity, K 
free surface, and ps is the stagnation pressure. 

In terms of our dimensionless variables, we rewrite (1) as 

is the curvature of the 

i(ir +v 2 )-K = 

(2) 

where 

, pb 2 

Z = Y 2 ( Ps-PB )• 

(3) 

The kinematic boundary conditions on the walls yield 


v = 0 on CD, 

(4 a) 

and 


i> = i/tana on AE. 

(4 b) 

Finally, we impose a contact angle fl at cj) = +1 by the relations 


v = ntan/i at 0 = 1, 

(5 a) 

v = — «tan(/l — a) at 0 = — 1. 

(5 b) 


In (5a) and (5b), the angle ft is assumed to be prescribed. 

This completes the formulation of the problem. We seek u — iv as an analytic function 
of / = cf> + ixp in the lower half plane ip < 0 satisfying equations (2), (4) and (5). 
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3 Numerical procedure 


In this section we derive a numerical scheme based on series truncation to solve the 
problem of Figure 1. The scheme supplements the work of Vanden-Broeck & Keller 
[12, 13], where solutions for a = n and a = n/2 were calculated by boundary integral 
equation methods. 

We first map the flow domain into the inside of a unit disk in the complex f-plane by 
the transformation 


/ = 


1 + 1 2 
2t 


( 6 ) 


This transformation maps the walls on the real diameter and the free surface on the 
unit circle (see Figure 3). The point B in Figures 2 and 3 is the point on the free surface 
where (/> = 0. We shall seek solutions which are symmetric with respect to B. 

We denote the complex velocity by £ = u — iv. We consider £ as a function of t. Since 
there are angled corners at t = +1, £ has singularities at these points. The appropriate 
singularities are 




c~(f+ 1) 


as 


+ 1 . 


(7) 


Furthermore as t —► 0, 



We represent the complex potential £ by the expansion 


( 8 ) 


C = (1 


* 2 \ 2 - 


21 

K 



00 

^ a n t n ~ l 


(9) 


The kinematic conditions (4) on the walls AE and CD imply that the coefficients a„ are 
real. It can be checked that (9) satisfies (4) and (5). Therefore, we can expect the series in 
(9) to converge for |f| < 1. The coefficients a„ must be determined to satisfy the condition 
(2) on the free surface ^4C. 

We use the notation t =\ t \ e m so that points on AC are given by t = e w , where 
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0 ^ o ^ 7t. Using (6), we find that 0 = 

dx . 0y 
df ~ *00 

we can then write the curvature as 


cos(o-) on the free surface. Since 
1 u + iv 

u — iv u 2 + v 2 ’ 


K _ 1 uv a —vu a 

sin(cr) f u 2 + v 2 ’ 

and Bernoulli’s equation (2) yields 


( 10 ) 


( 11 ) 


i(u 2 


2 , , 1 UV a — VU„ 

' v ) + _T— T“i / , , = £ 

sin(cr) ^u 2 +1 > 2 


( 12 ) 


on the free surface AC. 

We solve the problem approximately by truncating the infinite series in (9) after N — 2 
terms. Thus, we write 


£ ~ (1 — f 2 ) 2 " J2a n t n (13) 

' ' n=l 

and we introduce the N — 2 collocation points 

I = — 2. (14) 

Substituting f = e" 7 ' into (13) and taking the real and imaginary parts, we obtain the 
values of u and v at the mesh points (14) in terms of the coefficients a„. Similarly, we 
substitute t = e w into (13), differentiate with respect to a, and take the real and imaginary 
parts to obtain u„ and v a at the mesh points (14) in terms of a n . Substituting these values 
in (12) evaluated at the mesh points (14), we have N — 2 equations for the N— 1 unknowns 
i; and a n , n = 1, — 2. We obtain one more equation by fixing an extra parameter 5. 

Here we choose 

^a/(27u—a) 

<5 =( 15 ) 

a j 

The last equation is then given by (15), where the left-hand side is given. The particular 
choice (15) is motivated by the fact that (15) reduces to the parameters y used in Vanden- 
Broeck & Keller [12, 13] when a = n and a = n/2. This facilitates comparisons. For 
given values of a, /) and S, the system of N — 1 nonlinear algebraic equations with TV — 1 
unknowns is solved by Newton’s method. Once a solution was obtained for given values 
of a, /? and <5, a continuation method was used (i.e. the solution was used as an initial 
guess to compute a new solution for perturbed values of the parameters a, fl and 5 and 
so on). 


4 Discussion of results 

We used the numerical scheme described in § 3 to compute solutions for given values 
of a, fi and fi. The coefficients a„ were found to decrease very rapidly. For example, 
for a = /) = n/2 and <5 = —1.3, a\ « 0.93, a\\ » —0.36, « 2 i » 0.028, cqi » 0.000095 
and «65 ~ 0.000000072. Most of the results presented in this paper were obtained with 
N = 160. 
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Figure 4. Free surface profiles for (a) a = k/2 and p = n/2 with <5 = 1, (b) S = 0 and 

(c) <5 = S 0 = -1.3. 


Typical free surface profiles are shown in Figures 4-6. As 8 —► oo, the bubble tends 
to a circular arc. It is easily shown that it is consistent with (2) and (15). As 5 — ► oo, 
l, — > —oo and (2) implies that K —* Therefore, the bubble profile approaches an 
arc of a circle of radius — 1/£ as 8 —* oo. As <5 decreases the bubble elongates in the 
direction which bisects the angle a. When <5 reaches a certain value 8o(p, «), opposite sides 
of the bubble touch each other. For jS < a/2 this point of contact is at the intersection 
of the two walls (see Figure 6), while for ft > a/2 it is off the walls (see Figures 4 
and 5). 

The solutions for a = n were previously calculated elsewhere [3-6]. For p 4= n/2, they 
model the deformation of a bubble attached to a wall. For [} = n/2 they also represent 
half of free bubble. 

In Vanden-Broeck & Keller [13], solutions for a = n/2 and 8 > —1.24 were obtained 
by a boundary integral equation method. Here these results are extended up to the value 
8 = —1.3 at which opposite sides of the bubble touch at one point (see Figure 4). 

Finally, let us mention that our scheme can also be used to compute solutions for 
a > ji. It corresponds to a flow around a corner with a bubble attached at its apex. A 
typical profile is shown in Figure 7. We note that the velocity is finite everywhere along 
the streamline \p = 0. This is to be contrasted with the flow in the absence of a bubble, 
for which the velocity is infinite at the apex of the corner. 
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5 Conclusions 

We have used a series truncation method to calculate the deformation of a bubble in a 
corner flow. Our results include as particular cases those Vanden-Broeck & Keller [12, 13]. 
We were able to extend the calculations up to the limiting configurations, where opposite 
sides of the bubbles touch each other. 
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